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Abstract 

Background: In conditions of nitrogen limitation, Saccharomyces cerevisiae strains differ in their fermentation 
capacities, due to differences in their nitrogen requirements. The mechanisms ensuring the maintenance of 
glycolytic flux in these conditions are unknown. We investigated the genetic basis of these differences, by 
studying quantitative trait loci (QTL) in a population of 133 individuals from the F2 segregant population 
generated from a cross between two strains with different nitrogen requirements for efficient fermentation. 

Results: By comparing two bulks of segregants with low and high nitrogen requirements, we detected four 
regions making a quantitative contribution to these traits. We identified four polymorphic genes, in three of these 
four regions, for which involvement in the phenotype was validated by hemizygote comparison. The functions of 
the four validated genes, GCN1, MD53, ARG81 and BI03, relate to key roles in nitrogen metabolism and signaling, 
helping to maintain fermentation performance. 

Conclusions: This study reveals that differences in nitrogen requirement between yeast strains results from a 
complex allelic combination. The identification of three genes involved in sensing and signaling nitrogen and 
specially one from the TOR pathway as affecting nitrogen requirements suggests a role for this pathway in 
regulating the fermentation rate in starvation through unknown mechanisms linking nitrogen signaling to 
glycolytic flux. 
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Background 

Yeast strains are the main microorganisms used in fer- 
mentation process. During wine fermentation, yeast and 
principally Saccharomyces cerevisiae, consumes the sugars 
found in the grapes musts and converts them into alcohol, 
carbon dioxide and secondary-ends products that contrib- 
ute of wine character. To support yeast growth and enable 
it to perform these complex biochemical transformations 
a number of nutrients must be found in musts. The assim- 
ilable nitrogen is a key nutrient in the control of alcoholic 
fermentation and is consumed at the beginning of the 
process. Low nitrogen levels in musts may cause slow or 
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stuck fermentation [1]. During fermentation, the ethanol 
concentration increases and can destabilize cell mem- 
branes in a manner that may result in an inability to take- 
up nitrogenous compounds from the must [2]. Alcoholic 
fermentation occurs principally in these stressed condi- 
tions, so the ability of the yeast to maintain high levels of 
fermentation activity in such conditions is crucial to the 
outcome of the alcoholic fermentation. Wine yeast strains 
differ in their capacities to carry out fermentation in con- 
ditions of nitrogen limitation. These differences have been 
described as reflecting differences in the nitrogen require- 
ments of wine yeasts [3-5]. During alcoholic fermentation, 
such differences become visible after the yeast cells enter 
stationary phase in many cases though not all due to ni- 
trogen starvation (i.e. in a complete exhaustion of assimil- 
able nitrogen in must), through differences in the abilities 
of different strains to maintain a fermentation flux. Little 
is known about the mechanisms involved in controlling 
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fermentation rate under nitrogen starvation and the limit- 
ing factors involved. It has been suggested that control 
over sugar transport capacity affects the rate of fermenta- 
tion in nitrogen-starved cells. It has been shown that hex- 
ose carriers are targeted for degradation by endocytosis in 
such conditions, suggesting that carrier stability may affect 
the fermentation rate [6]. Yeasts are also subjected to the 
inhibitory effects of ethanol, which also decreases the fer- 
mentation rate [7]. Several studies have reported the 
genome-wide transcriptional response of yeast to nitrogen 
starvation in alcoholic fermentation conditions [4,5,8-12]. 
In a previous study, we showed that "low-nitrogen re- 
quirement strains" (LNR) strongly expressed biosynthetic 
genes, whereas "high-nitrogen requirement" (HNR) strains 
displayed a specific gene expression pattern, with the 
overexpression of stress genes. HNR strains seem to be 
more sensitive to nitrogen starvation, resulting in a 
stronger stress response and a lower fermentation rate. 
These differences in response are reminiscent of the dif- 
ferences in fermentation capacity described between 
sake and laboratory yeasts, resulting from differences in 
Riml5p function [13]. 

RimlSp is involved in nitrogen signaling downstream 
from TOR, which senses the nitrogen status of the cell 
and adapts cell metabolism to nutrient availability [14,15]. 
These previous studies provided the first demonstration 
that nitrogen signaling could affect fermentation rate. 
Similar mechanisms may contribute to differences in the 
nitrogen requirements of wine yeasts and would be con- 
sistent with differences in transcriptional patterns. How- 
ever, although physiological approaches have provided 
novel and relevant insight into the mechanisms associated 
with nitrogen requirements, the genetic variants under- 
lying phenotypic differences have yet to be identified. The 
use of QTL analysis to identify the genes underlying this 
variation should improve our understanding of the mech- 
anisms involved. 

Nitrogen requirement is a quantitative trait. It is there- 
fore possible to use a QTL approach to investigate the 
molecular basis of variation for this trait. QTL approaches 
have led to the detection of many genes, typically unlinked 
to each other [16]. In Saccharomyces cerevisiae, several 
QTL analyses have already been carried out, to identify 
genes involved in growth at high temperature [17,18], 
sporulation [19-22], cell morphology [23], drug sensitivity 
[24], ethanol tolerance and growth [22,23,25-27], floccula- 
tion [28], wine aroma production [29], amino acids con- 
sumption [30], and to decipher regulatory network 
variations [31,32]. These studies have shown some pheno- 
types to be highly complex. Some phenotypic variation 
may be accounted for by a single QTL with a major effect 
[33-36], or by several QTLs with minor effects [18,37,38]. 
It is thought that a large proportion of phenotypic vari- 
ation is accounted for by QTL with smaller effect sizes 



[39-41]. The complexity of the phenotype is accounted 
for not only by the number of genes underlying the vari- 
ation, but also by genetic interactions and epistasis be- 
tween loci [42]. 

We investigated the genetic basis of the variability of ni- 
trogen requirement, by establishing a genetic device based 
on bulk segregant analyses and using it to identify QTL. 
We identified four genes for which allelic variations be- 
tween parental strains were associated with differences in 
fermentative activity in a medium in which nitrogen was 
limiting. Interestingly, three of these four genes were 
found to be involved in nitrogen sensing and signaling. 

Results 

Screening of parental strains and constitution of the 
study population 

We investigated the genetic basis of variations of nitro- 
gen requirements in wine yeast, in two Saccharomyces 
cerevisiae enological strains, MTF2029 and MTF1782, 
characterized in a previous study as displaying extreme 
differences in fermentation ability in musts in which 
nitrogen was limiting [5]. Constant fermentation rate 
(CFR) determinations [43] indicated that strain MTF2029 
had a low nitrogen requirement (0.94 mg N g _1 C0 2 ), 
whereas MTF1782 had a high nitrogen requirement 
(2.5 mg N g _1 CO2). The nitrogen requirements of these 
strains were directly correlated with their fermentation 
capacity in a nitrogen-deficient medium (Figure 1). At the 
start of the stationary phase, when cell growth had 
stopped, fermentation rate declined differently in these 
two strains. Strain MTF2029 maintained a high fermenta- 
tion rate throughout stationary phase, whereas MTF1782 
displayed a large drop in fermentation flux, resulting in a 
longer fermentation time. The nature of the nitrogen 
source contained in the synthetic medium must be taken 
in consideration because it can impact the fermentation 
efficiency. Moreover, the fermentation rate is positively 
correlated with both the total amount of assimilable nitro- 
gen and the nitrogen uptake rate [12,31]. Indeed, yeast 
can display differences in nitrogen compounds consump- 
tion [30,44]. To avoid differences phenotypic caused by ni- 
trogen source utilization, we ensured that the two 
parental strains show the same consumption profile for 
nitrogen sources. These differences in fermentation rate 
were not associated with differences in the use of available 
nitrogen sources, as both strains used all the assimilable 
nitrogen during the growth phase. These two strains were 
therefore considered relevant for use in studies of the gen- 
etic basis of variation in fermentation capacity in condi- 
tions of nitrogen deficiency. The choice of these strains 
was also based on their ability to sporulate (approximate 
of 30% for the two strains), because this characteristic is 
required for the construction of a recombined population. 
We have estimated a spore viability of 33% for MTF2029 
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Figure 1 Comparison between the two parental strains (MTF2029 and MTF1782). Fermentation performances in nitrogen-deficient medium 
(SM100) (A). Comparison of the nitrogen requirements during fermentation. Nitrogen requirement determined by the constant fermentation rate 
method (B). 



strain and 44% for MTF1782 strain. A hybrid strain was ob- 
tained by crossing two parental haploid clones, 2029-C5 
and 1782-B1, in which the HO gene had been inactivated. 
The hybrid strain had a fermentation profile intermediate 
between those of the two parental strains (Additional file 1: 
Figure SI). The zygote was allowed to sporulate and the Fl 
haploid clones obtained were used for a second round of 
crosses, generating 133 F2 haploid segregants. 

Phenotyping of the segregant population 

We phenotyped the segregant population in nitrogen- 
limited fermentation conditions. We first determined the 



amount of CO2 produced after 89 hours of fermentation 
on a medium containing 120 mg of available nitrogen. 
We had previously checked that this measurement was 
representative of the fermentative capacity of the strain 
and correlated with nitrogen requirement determina- 
tions by the CFR procedure. The whole population was 
phenotyped in these conditions and the results are 
shown in Figure 2. This phenotypic trait displays a low 
transgression level with only 8 segregants in negative 
transgression corresponding at 6% of transgressive segre- 
gation in the progeny. More the control of variability re- 
vealed a heritability of 98% for the character indicating 




F2 segregants population C0 2 to 89 h fermentation (g L J } 



Figure 2 Distribution of the amount of C0 2 released at 89 hours in the segregating population in nitrogen-deficient medium (SMI 20). 

Fermentation rate was measured as the amount of C0 2 released at 28°C. Segregants were sorted according to cumulative C0 2 release (A) and 
the fermentation rate histogram (B). For both representation, purple and orange code represents the position of 1782-B1 and 2029-C5 parental 
strain, respectively. 
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that phenotypic variation resulted from a strong genetic 
control. The continuous distribution of "fermentation 
capacity" in the population, suggests that this character 
is quantitative in nature. Phenotype distribution was also 
consistent with the phenotype being determined by sev- 
eral genes. Comparison of the parental strains with the 
segregant population revealed that only one segregant 
had a greater fermentation capacity than the parental 
2029-C5 strain, whereas 16 segregants had a lower fermen- 
tative capacity than the 1782-B1 strain. Thus, the allelic 
combination responsible for the low-nitrogen requirement 
phenotype is rare. Conversely, the fermentative capacity of 
the segregants was lower than that of strain MTF1782, indi- 
cating that there are various negative alleles in the genomes 
of these two parental strains that can be combined. 

We generated a large population of F2 segregants 
(133), for the detection of combinations of several QTL 
explaining this complex phenotype. We overcame the 
difficulties associated with the genotyping of such a 
large population by adopting a bulk segregant analysis 
strategy. The most distant phenotypes were selected 
and combined in two pools of 15 segregants each. Each 
strain chosen for the two bulks the "low-nitrogen re- 
quirement" (LNR-bulk) and "high-nitrogen requirement" 
(HNR-bulk) bulks was characterized in greater detail by 
online fermentation monitoring. The segregants within 
each pool had similar fermentation profiles, whereas the 
two pools had contrasting fermentation profiles (Figure 3). 
The differences between the two bulks were most clearly 
visible during the onset of stationary phase (mid-fermen- 
tation at 50 hours of fermentation) and at the end of 
fermentation. 



Identification of the QTL region by bulk segregants 
analysis (BSA) 

For the detection of genomic regions involved in the 
"fermentation capacity" trait, we hybridized genomic 
DNA from the two pools of 15 segregants and the two 
parental strains to an Agilent 8xl5K custom isothermal 
array containing 6318 SNPs. The array was obtained 
from a comparison of the genomic sequences of P3-D5 
and RM11, strains genetically similar to the 2029-C5 
and 1782-B1 parental strains, respectively (see Methods). 
The hybridization signals for the 1900 significant oligo- 
nucleotides are plotted along the chromosomes for the 
two pools of segregants and the two parental strains 
(Figure 4A). There are large and significant (indicated by 
green spots) differences between the two bulks in several 
regions. For some regions, the maximum hybridization 
intensities indicate the almost exclusive presence of a 
single parental allele in one bulk. For example, this was 
observed once on the chromosome IV homologous to 
the S288C chromosome (1,163,169 to 1,214,335 bp), 
whereas the HNR-bulk displayed the same signal as the 
1782-B1 parent, and three times on the chromosome VII 
homologous to the S288C chromosome (98,215 bp to 
147,601 bp; 674,690 bp to 691,654 bp; 730,853 bp to 
765,691 bp) whereas the HNR-bulk displayed the same 
signal as the 2029-C5 parent for the first two peaks and 
the signal of the 1782-B1 parent for the third peak. The 
hybridization profiles suggest that 23 regions differed 
significantly between the two bulks. However, the small 
population used for each bulk may have resulted in local 
variations in the hybridization signal, leading to the de- 
tection of false-positive peaks. We reduced this risk to 




Time (h) 

Figure 3 Fermentation profiles obtained for strains from each bulk in nitrogen-deficient medium (SMI 20) at 28°C. For fermentation 
representation, purple and orange code represents the HNR-bulk and the LNR-bulk, respectively. 
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Figure 4 Quantitative allele frequency measurement in DNA pools. Genotyping of parental strains (2029-C5 in brown and 1782-B1 in deep 
purple), and two segregating pools with low (LNR orange) and high (HNR in purple) nitrogen requirements. The 2029-C5 allele enrichment of the 
pools is indicated by deviations above 0 and 1782-B1 allele enrichment is indicated by deviations below 0. Regions with significant differences in 
allelic frequencies between the two pools are indicated by light green dots. RM1 1 probes were positioned on the S288C genome by Blast. Four 
regions displaying large differences between the two pools were chosen for further analysis: Al, A2, B and C. 



obtain false positive by restricting arbitrarily our analysis 
to regions for which the differences between the 
hybridization signals obtained for the two bulks for each 
allelic probe were significant (as determined with a 20- 
SNP window f-test) and greater than half the difference 
between the parents. Given the high number of peaks, 
genomic regions exhibiting the highest differences be- 
tween each bulk are likely to have the highest impact on 
the phenotype. The resulting four QTL regions were 
chosen for analysis (Figures 4A and 4B): two on the 
chromosome homologous to S288C chromosome VII, 
which were approximately 17 and 50 kb long (regions 
Al and A2), one 24 kb-long region on the chromosome 
homologous to S288C chromosome XIII (region B) and 
one 14 kb-long region on the chromosome homologous 
to S288C chromosome XIV (region C). 

For regions Al, B and C, the bulk with the low nitrogen 
requirement (LNR-bulk) contained the alleles of the par- 
ent with the opposite pattern of nitrogen requirement, 



whereas the bulk with the high nitrogen requirement 
(HNR-bulk) contained specific markers of the 2029-C5 
parent. By contrast, the LNR-bulk and the HNR-bulk con- 
tained 2029-C5 and 1782-B1 alleles, respectively, in the 
A2 region. These results suggest that, for three of the four 
QTL regions identified, the allele increasing fermentation 
rate originates from the strain with the lowest fermenta- 
tion capacity (1782-B1), whereas only one allele from the 
strain with the highest fermentative capacity (2029-C5) 
had a positive effect on fermentation kinetic at one locus. 

Characterization of polymorphism for the candidate 
genes 

We examined the four QTL regions, to identify the 
best candidate genes (Additional file 2: Figures S2 and 
Additional file 3: Figure S3). We chose four genes involved 
in nitrogen metabolism: MDS3, GCN1, ARG81, and BI03. 
MDS3, which is located in QTL region A2, encodes a 
component of the TOR pathway, and GCN1, from the 
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same region, is a regulator Gcn2p kinase involved in gen- 
eral amino-acid control [45]. These two genes are involved 
in nitrogen sensing and signaling in response to nitrogen 
availability. ARG81, located in the QTL region B, is a rele- 
vant candidate gene because it encodes a transcription fac- 
tor that regulates gene expression in response to arginine 
availability. The BI03 gene is located in QTL region C 
and encodes a protein involved in biotin synthesis. 

We investigated whether polymorphism of the coding 
regions of these genes could account for phenotypic vari- 
ation, by sequencing the four genes of the two parental 
strains (Figure 5A). We detected several nonsynonymous 
allelic variants (Figure 5B) that could potentially account 
for differences in nitrogen requirement. As phenotypic 
differences may also reflect differences in gene expression, 
we also compared the expression levels of these genes, as 
determined in a previous study [5]. No significant differ- 
ence was found between MTF1782 and MTF2029 for the 
expression of these four genes, suggesting that the varia- 
tions in expression were not involved. 

Functional analysis 

For each candidate gene, we constructed two hemizygotic 
strains, by crossing one parent bearing an inactivated form 
of the gene with the other parent containing a functional 
form. Each hemizygote was phenotyped on nitrogen- 
deficient medium (SM100) and the fermentation kinetics 
of the strains were compared (Figures 6 and 7). 

MDS3 was the only gene tested for which a hemizy- 
gous construction with an allele from strain MTF2029 
(1782-BlAx2029-C5) resulted in a higher fermentation 
capacity than a hemizygous construction with an allele 
from MTF1782 (Figure 6A). The difference between the 
two hemizygous constructions corresponded to about 
7% of the difference between the two parental strains in 



stationary phase (70 hours into the fermentation) and 
15% of total fermentation time. MDS3 was also the only 
gene for which a hemizygous strain had a higher fer- 
mentation capacity than the hybrid. The simultaneous 
expression of the two alleles in the hybrid thus has a 
negative effect on fermentation, resulting in a lower fer- 
mentation capacity than the expression of the allele ori- 
ginating from MTF2029 alone. Studies of hemizygotes 
for the GCN1 gene showed that the allele from parent 
1782-B1 conferred a higher fermentation capacity than 
the allele from 2029-C5 (Figure 6B). Total fermentation 
time differed by about 20% between these two hemizy- 
gotes. The effects of the GCN1 and MDS3 genes were 
interesting, because both these genes are located in 
QTL region A2. Indeed, they are separated by 2364 
nucleotides but have opposite effects on fermentative 
capacity on nitrogen-deficient medium (Figure 6C). The 
determination of a possible complex genetic interact- 
ion between alleles MDS3 and GCN1 requires further 
characterization of this QTL structure. The reciprocal 
hemizygosity tests performed for each single gene did not 
permit to assess the combined effects of alleles. In the par- 
ental strain 2029-C5, the MDS3 gene has a positive impact 
on this phenotype, whereas the GCN1 gene has a negative 
effect (Figures 6A and 6B). For the parental strain 1782- 
Bl, phenotypic comparison and nucleotidic variation 
between different strains show that MDS3 is a recessive 
allelic form and GCN1 is a dominant allelic form. 

The differences in fermentation kinetics between the two 
hemizygous constructions (Figure 7A) confirmed the in- 
volvement of ARG81 in fermentation capacity. The allele 
from the parental strain 1782-B1 conferred a better fermen- 
tative performance than the allele from parental strain 
2029-C5, with fermentation time about 29% shorter for the 
1782-B1 allele. However, BI03 hemizygous constructions 
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Figure 5 Nucleotide sequence variation for the four candidate genes. Nucleotide sequence variation between the two parental strains 
(2029-C5 and 1782-B1) (A). Amino-acid differences between strains (B). 
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had the most significant effect on fermentation rate under 
nitrogen limitation conditions (Figure 7B). The difference 
between the two hemizygotes corresponded to about 20% 
of the difference between the parental strains for fermenta- 
tion rate during stationary phase and total fermentation 
time. The hemizygous phenotypes of ARG81 and BI03 
were consistent with the reverse distribution of markers in 
these two QTL regions: i.e. the HNR-bulk contains specific 
markers of 2029-C5. 



The total fermentation rate and the exponential 
growth phase were unaffected by the inactivation of the 
four genes, in all the constructions tested. A comparison 
of cell growth revealed no significant difference between 
the two hemizygous constructions for MDS3, BI03 and 
ARG81 (Figure 8). For the GCN1 gene, the growth pro- 
files of the hemizygotes showed that the allele from 
1782-B1 conferred higher levels of growth, resulting in a 
higher biomass, than the allele from 2029-C5. 
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These results confirm the involvement of these four 
genes in fermentative capacity in conditions of nitrogen 
limitation. The phenotypes obtained for the various 
hemizygotes remained very different from the pheno- 
types of the haploid parents. These findings highlight 
the complexity of this technological phenotype and dem- 
onstrate that the phenotype of the two haploid parents 
results from the allelic combinations of several genes, 
which may also interact. 

Discussion 

In this study, we constructed and characterized a recom- 
bined F2 segregant population and used a BSA approach 
to identify genes involved in the variation of fermenta- 
tion capacity under conditions of nitrogen limitation. 
This led to the detection of four genomic regions in- 
volved in "low" or "high" nitrogen requirement pheno- 
types. For three of these four regions, we were able to 
demonstrate that allelic variations of four genes had an 
impact on fermentation performance. Unexpectedly, for 
three of the four genes, the allele with a positive effect 
on fermentation rate originated from the strain with a 
low fermentation capacity in such conditions (1782-B1), 
rather from the strain with a high fermentation capacity, 
2029-C5. This finding of multiple genes with positive ef- 
fects in strains with a low fermentation capacity is 



consistent with the complex character of the trait, as 
suggested by the distribution of "nitrogen requirement" 
values in the recombined population. As we identified 
only one gene with a positive effect originating from 
2029-C5, other genes not detected in this first analysis 
are probably involved in this phenotype, as suggested by 
the numerous peaks in the hybridization profile, corre- 
sponding to putative QTL regions. Our ability to detect 
more QTL is, of course, limited by the size of the segre- 
gant population, which also imposed constraints on the 
number of segregants per pool that we could analyze 
and by the phenotypic variation between the two se- 
lected bulks. The phenotyping of a larger population 
might lead to the identification of other genes involved 
in this complex phenotype. 

One of the most interesting findings of this study is the 
organization of QTL region A2 into a complex architec- 
ture, with two alleles of MDS3 and GCN1 having opposite 
effects on the phenotype in the 2029-C5 strain. The 
MDS3 allele has a positive effect on fermentation rate 
whereas the GCN1 allele has a negative effect, and both 
genes are involved in mechanisms responding to nitrogen 
availability. The QTL was detected by hybridization as a 
region from 2029-C5 with a positive impact on fermenta- 
tion rate (enriched in the LNR pool). As this QTL con- 
tained the two genes, the impact of GCN1 allelic variation 
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would be expected to be much weaker than that of MDS3. 
However, this supposed difference in power between the 
two genes was not observed in functional analysis based 
on reciprocal hemizygote tests, in which both alleles had a 
moderate to weak effect, suggesting that this assay may 
underestimate the effect of the MDS3 allele. While the im- 
pact on the phenotype is low, MDS3 and GCN1 have an- 
tagonistic effects on the phenotype. Moreover these two 
genes have functional relationship. Physical and functional 
proximity of these genes suggest a possible interaction be- 
tween alleles. This would be consistent with a moderate 
effect of the two alleles on the phenotype. This combin- 
ation could provide a selective advantage to the strain and 
will justify the association of alleles with antagonistic ef- 
fects. We could also suppose that this phenotype is the re- 
sult of interaction with other genes located in other 
regions. More detailed characterization of the genes will 
be required to determine the impact of allelic variation. 
Three of the four genes identified as candidate genes in 
this study are directly involved in nitrogen metabolism, in 
sensing and signaling or regulating gene expression in re- 
sponse to nitrogen status. The remaining gene, BI03, is 
involved in biotin metabolism, which is not directly linked 
to nitrogen metabolism. 

These four genes harbor mutations predicted to result 
in amino-acid substitutions potentially affecting the ac- 
tivity of the encoded protein. All non-synonymous SNPs 
were analyzed by SIFT (Sorting Intolerant From Toler- 
ant). This program is a tool based on sequence hom- 
ology for predicting the probability that a mutation 
might be deleterious and affects the protein function 
[46]. Phylogenetic analysis (Additional file 4: Figure S4) 
showed that the two allelic forms of the different genes 
were relatively similar. Mds3p displays four nonsynon- 
ymous mutations with respect to the two parental strains, 
and two mutations are positioned in N-terminus protein. 
According to SIFT analysis, these mutations should not 
affect the protein function in a deleterious manner. Phylo- 
genetic analysis suggested that two of the variants present 
in the 2029-C5 parental strain were ancestral whereas the 
other two were ancestral in the 1782-B1 strain (Figure 5B). 
For Gcnlp, phylogenetic analysis showed that the two 
nonsynonymous mutations in 1782-B1 corresponded to 
the ancestral form of the allele (Figure 5B). The second 
mutation is particularly interesting because it affects an 
Armadillo -like helical domain of the protein critical for 
its spatial conformation [47]. In this case, SIFT analysis 
indicate that the two mutations may impact the protein 
function. Comparative sequence analysis of the Arg81p 
protein showed that the two parents each have one mu- 
tation with respect to the ancestral allelic form (Figure 5B). 
The Bio3p protein displays one nonsynonymous substi- 
tution, for which it is difficult to confirm the ancestral 
allelic form. This substitution is located in a pyridoxal 



phosphate-dependent transferase domain of the protein. 
For ARG81 and BI03, SIFT analysis did not reveal pos- 
sible deleterious mutation. 

For the four genes studied, the impact of the mutations 
identified on the activity of the proteins was unknown and 
difficult to infer from our data. MDS3 is involved in the 
TOR signaling pathway and was shown to function as a 
positive regulator acting on TAP42 [48]. TOR senses the 
nutrient status of the cell (specifically amino acids) and 
coordinates the cellular response through the control of 
protein synthesis and ribosome biogenesis, which are 
stimulated when TOR is active, whereas growth stops and 
stress response are triggered by nitrogen starvation and 
TOR inactivation. A previous transcriptomic comparison 
of the parental strains revealed differences in gene expres- 
sion, such as higher levels of ribosomal gene expression in 
MTF2029 and of stress gene expression and protein deg- 
radation in MTF1782, consistent with a higher level of 
TOR pathway activity in the MTF2029 strain compared to 
MTF1782 [5]. MDS3 expression did not differ significandy 
between the two strains but, as MDS3 activates the TOR 
pathway, the variations of gene expression observed are 
consistent with higher activity of Mds3p in MTF2029. 
Moreover, this result is reminiscent of the impact of vari- 
ation of the RIM1S gene, encoding a TOR-controlled PAS 
kinase, on the fermentation rate of sake yeasts [13]. This 
previous study showed that changes in RIM1S function 
prevented the entry of the cells into quiescence during 
starvation and led to the maintenance of high rates of gly- 
colysis. Indeed, higher levels of MDS3 activity could 
modulate TOR activity but the mechanisms by which such 
an increase in TOR activity could lead to an increase in 
fermentation rate are unclear. Studies suggest that TOR 
can impact glycolytic flux by a modulation of the PKA ac- 
tivity through different intermediates such a Gcn4p [49], 
FGM pathway [50] or ammonium permeases which act as 
sensor [51]. 

GCN1 is known to be a positive regulator of general 
amino-acid control (GAAC) [45]. It is required for the 
activity of GCN2, which triggers the preferential transla- 
tion of mRNA encoding the GCN4 transcription factor, 
which controls the expression of genes encoding pro- 
teins involved in amino-acid biosynthesis [45]. Gcnlp 
forms a complex with Gcn20p. This complex is not re- 
quired for GCN2 activation, but increases the phosphor- 
ylation of eIF2a by Gcn2p [45]. However, GCN2 regulation 
is not thought to be functional in nitrogen-starved cells 
[52]. Two mutations of the GCN1 gene can be used to 
distinguish between the two parental allelic forms. One 
of these mutations affects an Armadillo repeat domain 
involved in the spatial conformation of this protein. It 
remains unclear how variations affecting Gcnlp causes 
changes in fermentation rate. It is thought that changes 
in the conformation of the protein may influence the 
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formation of the Gcnl/Gcn20 complex and GCN2 acti- 
vation. A perturbation of GNC2 activation could impact 
the activation of GCN4 and intracellular amino-acid 
pools. ARG81 is a transcription factor included in a pro- 
tein complex (Arg81p, Arg80p, Mcmlp, Arg82p), control- 
ling the expression of genes involved in the anabolism and 
catabolism of arginine [53]. Indeed, mutations in ARG81 
gene may lead to perturb the arginine metabolism and 
change the cellular pool of nitrogenous compounds [54]. 
It is not possible to infer the impact of the differences be- 
tween the two ARG81 forms on protein activity. However, 
it is possible that an increase in ARG81 activity might lead 
to changes to the cellular amino-acid pool, with an impact 
on protein synthesis or nitrogen signaling. 

BI03 encodes a DAPA aminotransferase involved in 
biotin biosynthesis [55]. The BI03 gene is not directly 
connected to nitrogen metabolism. However, biotin is 
involved in carboxylation reactions, some of which are 
involved in nitrogen metabolism, such as the reactions 
catalyzed by urea carboxylase or pyruvate carboxylase 
(generating oxaloacetate, a precursor of a-ketoglutarate 
and aspartic acid). Changes in biotin availability might 
therefore have an effect on the equilibrium of the nitro- 
gen pool. Alternatively, its is notable that BI03 uses 
SAM (S-Adenosyl methionine) as substrate containing a 
nitrogen group and its activity might impact on SAM 
pool which is involved in other cellular metabolisms. 
The only mutation found in this gene affected the pyri- 
doxal phosphate-dependent transferase domain. This do- 
main is involved in key functions, such as the transfer of 
nitrogenous groups. Mutations affecting this domain 
may affect the concentration of biotin, which is required 
by the cell. 

Conclusion 

In conclusion, we show here that fermentation perfor- 
mances in conditions of nitrogen limitation result from 
the cumulative effects of multiple alleles, suggesting that 
this phenotype is highly complex. This QTL study led to 
the identification of three candidate genes with functions 
associated with the regulation of nitrogen metabolisms, 
nitrogen sensing or signaling: GCN1, ARG81 and MDS3. 
These findings highlight the role of nitrogen signaling in 
the control of glycolytic flux in nitrogen starvation and 
support the hypothesis that the TOR pathway plays a key 
role in controlling fermentation capacity in nitrogen- 
starved cells, consistent with previous observations [56]. 
These finding are consistent with the previous data indi- 
cating that the two strains did not have differences in pro- 
tein synthesis capacity in such starved conditions. 
Additional studies are required to identify the precise 
mechanisms by which variations of TOR signaling modu- 
late fermentation flux. 



Methods 

Construction of the parental strains 

We used two enological Saccharomyces cerevisiae strains 
in this study: MTF2029 and MTF1782. The nitrogen re- 
quirements and phenotypic characteristics of these two 
strains were determined in a previous study [5]. Both 
these strains are homothallic HO/HO diploids. The HO 
gene was disrupted in both strains, to obtain haploid 
clones for the crossing experiments. Yeasts were trans- 
formed, as described by Guldener et al. [57], with pUG6 
(KANMX6 cassette) for MTF2029 and pAG25 (NATMX4 
cassette) for MTF1782. For the disruption of one copy of 
the HO gene, we used the same primers for both parents. 
The two custom-made 60-mer primers used began with 
40 nucleotides identical to the upstream or downstream 
region of the HO gene, followed by 20 nucleotides for 
amplification of the KANMX6 and NATMX4 cassettes 
(Forward primer: ATGCTTTCTGAAAACACGACTATT 
CTGATGGCTAACGGTGCTTCGTACGCTGCAGGTC 
and Reverse primer: TTAGCAGATGCGCGCACCTGCG 
TTGTTACC ACAACTCTTTAGTGGATCTGATATCAC 
CTA). Sporulation was allowed to occur and we then dis- 
sected transformants displaying integration of the drug 
resistance cassette. We collected haploid ho spores, which 
we grew on YPD (yeast extract, peptone and dextrose) 
plates supplemented with the corresponding antibiotic 
(G418 at 200 ug/ml for MTF2029 and cloNAT at 
200 (rg/ml for MTF1782J. Two haploid ho spores were 
selected: 2029-C5 from MTF2029 and 1782-B1 from 
MTF1782. A zygote was created by crossing 2029-C5 
(Mat a) and 1782-B1 (Mat a) and tested on YPD plates 
supplemented with G418 and cloNAT. 

Construction of the segregant population 

We generated an F2 population of segregants, as a 
means of obtaining more recombinants. The Fl segre- 
gant population was obtained by tetrad dissection of the 
parental cross 2029-C5xl782-Bl. An F2 diploid popula- 
tion was then constructed by randomly crossing Fl seg- 
regants from different asci. The F2 diploids were then 
allowed to sporulate and the spores were collected by 
tetrad dissection. The ploidy status of the F2 segregants 
was checked in a mating test [26]. 

Fermentation conditions 

Phenotypic parameters were measured in batch fermen- 
tations in a synthetic medium (SM) as described by Bely 
et al. [58], mimicking a natural must contain: glucose 
(200 g liter 1 ), malic acid (6 g liter 1 ), citric acid (6 g liter -1 ), 
MgS0 4 (250 mg liter" 1 ), KH 2 P0 4 (750 mg liter" 1 ), CaCl 2 
(155 mg liter" 1 ), NaCl (200 g liter" 1 ), K 2 S0 4 (0,5 g liter" 1 ), 
vitamins and oligoelements mixtures. SM100 and SM120 
was supplemented respectively with a final concentration of 
100 and 120 mg liter" 1 assimilable nitrogen corresponding 
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to ammonium salt and a mixture of 19 amino acids (L-pro- 
line, L-glutamine, L-arginine, L-tryptophan, L-alanine, 
L-glutamic, L-serine, L-threonine, L-leucine, L-aspartic 
acid, L-valine, L-phenylalanine, L-isoleucine, L-histidine, 
L-methionine, L-tyrosine, L-glycine, L-lysine and L- 
cysteine). Batch fermentations were performed in a 1.2- 
liter fermenter and microfermenter (300 ml), with mixing 
by a magnetic stirrer (500 rpm), and airlocks to maintain 
anaerobiosis. We assessed the fermentation performances 
of the parental strains and hemizygotes in SM100 at 24°C, 
for precultures and fermentations. We carried out segre- 
gant analysis at 28°C in SM120, to make it possible to 
identify differences in fermentation performance more 
rapidly. 

Measurement of phenotypic parameters 

Two phenotypic measurements based on the CO2 re- 
leased during fermentation were obtained. The first 
phenotypic measurement was obtained at the start of 
the stationary phase (89 hours into the fermentation on 
SM120). Microfermentations were not monitored and 
fermenter mass was assessed manually. This measure- 
ment made it possible to distinguish rapidly between 
strains on the basis of their fermentation capacity on 
entry into stationary phase. This parameter was strongly 
correlated with the fermentation performance of the 
strains on nitrogen-deficient medium. 

We studied the kinetics of fermentation by online 
monitoring (fermentation batch in 1.2 liter). The amount 
of CO2 released during fermentation was calculated from 
automatic measurements (taken every 20 min) of fer- 
menter mass [59]. This fermentation monitoring method 
was validated in a previous study [60]. The rate of CO2 
production was calculated by polynomial smoothing of 
the last 10 measurements of fermenter weight loss. The 
many acquisitions of data for the mass and the precision 
of weighing (0.1 to 0.01 g) made it possible to calculate 
the rate of CO2 production with a high level of precision 
[58]. The high frequency of online measurements of CO2 
production (one measurement every 20 minutes) made 
it possible to calculate the rates of CO2 production by 
sliding-window second-order polynomial fitting in a 
custom-developed Labview application. This measure- 
ment made it possible to determine the overall rate of 
fermentation during stationary phase and to compare 
segregants at the same stage of fermentation. 

Reciprocal hemizygosity analysis 

This technique was used for the identification of allele- 
specific contributions to the phenotype [61]. We used 
the same primers for each candidate gene and the same 
transformation strategy as for the two parental strains. 
We inactivated one copy of the candidate gene with the 
hygromycin selection cassette. We then carried out a 



specific PCR to check that the gene was indeed inacti- 
vated. The primers used for cassette integration are indi- 
cated in the supplementary data (Additional file 5: 
Figure S5). The contribution of the alleles to the pheno- 
type was analyzed by comparing fermentation kinetics in 
SM100 (1.2-liter fermenters). 

Microarray design 

For the detection of genetic variation in the segregant 
population, we designed a DNA microarray with 
isothermal-melting probes, as described by Gresham 
et al. [62]. The genome sequences of Saccharomyces cer- 
evisiae RM11 and Saccharomyces cerevisiae P3-D5 (re- 
lated to MTF1782 and MTF2029, respectively) were 
used to construct a set of 6,318 pairs of isothermal 
probes. Biallelic positions at least 100 bp apart were 
chosen. We also included 300 replicates and 1000 invari- 
ant control probes in the array. The primers were de- 
signed with primer3 [63], to ensure hybridization at 50°C. 
The array design is available on GEO under the accession 
number GPL18217. 

DNA extraction, labeling and hybridization conditions 

For bulk segregant analysis, we selected two pools of 15 
segregants each. The first pool contained the 15 strains 
with the best fermentation performances in nitrogen de- 
ficiency conditions and the second pool contained the 
15 segregants with the poorest fermentation perform- 
ance in such conditions. Each pool was treated separ- 
ately. For each segregant, a culture in 30 ml of YEPD at 
28°C was prepared and the number of yeast cells present 
in the culture was determined by counting with an elec- 
tronic particle counter (Multisizer 3 counter; Beckman 
Coulter). For each pool, we isolated genomic DNA from 
all the segregants, mixed in equal proportions (5 x 10 9 
cells for each segregant), on Qiagen Genomic-tip 100/G 
columns. We also prepared genomic DNA from the par- 
ental strains, beginning with the same number of cells. 
For labeling and hybridization conditions, we used a 
modified version of the protocol described by Gresham 
et al. [62]. Genomic DNA was digested, fragmented by 
sonication and purified with a Qiagen Purification PCR 
kit. We labeled 1 ug of fragmented genomic DNA with 
Cy3 (for the two parental strains and the two bulks) or 
Cy5 (for a mixture of the DNA from the two parental 
strains) and carried out Agilent oligonucleotide array-based 
CGH for genomic DNA analysis, in accordance with the 
manufacturer's instructions (ref G4410-90010). Microarrays 
were hybridized at 59°C for 17 hours, then washed by 
immersion in two chip wash buffers: a low-stringency 
buffer (CGHmix2), followed by a high-stringency buffer 
(CGHmixl). Microarrays were scanned with a Genepix 
4000B scanner (Axon Instruments Inc.). 
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Statistical analysis 

Statistical analysis was performed with R software, version 
2.15.2 [64]. The phenotype heritability H 2 was calculated 
as previously described [65], i.e. H 2 = ((Var seg - Var env )/ 
Var seg ) xlOO, where Var env is the pooled variance among 
parental measurement and Var seg is the variance among 
phenotype values for the segregants. Transgressive segre- 
gation was defined as in [65,66] by the number of segre- 
gants whose phenotype level lay at least 2a higher than 
the mean phenotype level of the higher parent or 2a lower 
than the mean phenotype level of the lower parent; o is 
the pooled standard deviation of the parents. 

The Agilent 8xl5K array was imported into R software 
with the limma package [67]. For each probe of each 
block, the log2 (green signal/red signal) ratio was calcu- 
lated and centered. The red signal (Cy5) corresponded 
to the reference signal and the green signal (Cy3) corre- 
sponded to the signal for each bulk and each parent. Dif- 
ferences between log ratios corresponding to the P3-D5 
and RM11 alleles at biallelic loci were calculated and 
normalized, such that the parental strain MTF2029 had 
a mean value of 1.5 and the parental strain MTF1782 
had a mean value of -1.5. This criterion was used to 
analyze the differences between bulks. A first selection 
of probes was performed by carrying out a one-tailed t 
test based on parental strain criteria. Approximately 
1900 specific probes were selected after Benjamini- 
Hochberg correction for multiple testing [68] (adjusted 
j?-values < 0.05). We tried to identify allelic positions dif- 
ferentiating the HNR-bulk from the LNR-bulk for these 
1900 specific probes, by carrying out t tests on sliding 
sliding windows (20 probes per window) along the 
length of the genome. Benjamini-Hochberg correction of 
the ^-values was performed. Only probes with adjusted 
^-values < 0.01 at positions at which differences between 
bulks were greater than one third the differences be- 
tween the parents were retained (328 probes investi- 
gated). The complete array dataset (raw and processed 
data) is available from the Gene Expression Omnibus 
database under accession number GSE54389. 

Sequence analysis 

For each candidate gene, we compared gene and protein 
sequences, to identify nonsynonymous mutations distin- 
guishing between the two parental alleles. The poly- 
morphic change was analyzed by SIFT (Sorting Intolerant 
From Tolerant) [46]. For each allelic form, we ana- 
lyzed its distribution in the genomes of strains avail- 
able from the Saccharomyces Genome Database (SGD, 
http://www.yeastgenome.org) and strains from the Sac- 
charomyces Genome Resequencing Project (SGRP, 42 
strains) [69,70]. Phylogenies were inferred with MEGA 
[71], by the maximum likelihood method, based on 
the Kimura two-parameter model [72]. The trees with 



the highest log likelihood are shown. The trees are 
drawn to scale, with branch lengths proportional to the 
number of substitutions per site. 
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Additional file 2: Figure S2. QTL region map for the region Al (A) and 
the region A2 (B). For the region A2, the map is divided into two parts. 

Additional file 3: Figure S3. QTL region map, for the region B (A) and 
the region C (B). 

Additional file 4: Figure S4. Molecular phylogenetic tree for the four 
candidate genes (MDS3: A, GCN1: B, ARG81: C, BI03D). Evolutionary 
history was inferred by the maximum likelihood method, based on the 
Kimura 2-parameter model and using 43 nucleotide sequences from the 
available genome sequences [69] (SGRP2), 70 (SGRP1), [72]. 

Additional file 5: Figure S5. Primers used for cassette integration for 
gene inactivation in the parental strains. 



Competing interests 

The authors declare that they have no competing interest. 
Authors' contribution 

BB and JLL conceived designed and coordinated the study. CB performed 
experiments, QTL analysis, candidate gene search, and molecular 
constructions for genes validation. FB and JLL designed the isothermal 
genotyping array. IS performed statistical analysis (QTL). CB, JLL and BB wrote 
the manuscript. All authors read and approved the final manuscript. 

Acknowledgments 

This work was supported by the ANR project ALIA 2009. We thank Martine 
Pradal and Pierre Delobel for technical assistance. 

Author details 

1 1NRA, UMR1083 Science pour I'CEnologie, 2 Place Viala, F-34060 Montpellier, 
France. 2 Montpellier SupAgro, UMR1083 Science pour I'CEnologie, 2 Place 
Viala, F-34060 Montpellier, France. 3 Universite Montpellier 1, UMR1083 
Science pour I'CEnologie, 2 Place Viala, F-34060 Montpellier, France. 

Received: 31 January 2014 Accepted: 10 June 2014 
Published: 19 June 2014 

References 

1. Blateyron L, Sablayrolles JM: Stuck and slow fermentations in enology: 
statistical study of causes and effectiveness of combined additions of 
oxygen and diammonium phosphate. J Blosci Bioeng 2001, 91:184-189. 

2. Bauer FF, Pretorius IS: Yeast stress response and fermentation efficiency: 
how to survive the making of wine - a review. 5 Afr J Enol Vitic 2000, 
21:2751. 

3. Bely M, Sablayrolles JM, Barre P: Description of alcoholic fermentation 
kinetics: its variability and significance. Am J Enol Vitic 1990, 159:25-32. 

4. Gutierrez A, Chiva R, Sancho M, Beltran G, Arroyo-Lopez FN, Guillamon JM: 
Nitrogen requirements of commercial wine yeast strains during 
fermentation of a synthetic grape must. Food Microbiol 2012, 31(1):25-32. 

5. Brice C, Sanchez I, Tesniere C, Blondin B: Assessing the mechanisms 
responsible for differences in nitrogen requirements between 
Saccharomyces cerevisiae wine yeasts in alcoholic fermentation. 
Appl Environ Microbiol doi:10.1 128/AEM.03856-13. 

6. Salmon JM: Effect of sugar transport inactivation in Saccharomyces 
cerevisiae on sluggish and stuck enological fermentations. 

Appl Environ Microbiol 1989, 55(4)353-958. 

7. Ansanay-Galeote V, Blondin B, Dequin S, Sablayrolles JM: Stress effect of 
ethanol on fermentation kinetics by stationary-phase cells of 
Saccharomyces cerevisiae. Biotechnol Lett 2001, 23(9):677-681. 



Brice et al. BMC Genomics 2014, 15:495 
http://www.biomedcentral.com/1471-2164/15/495 



8. Rossignol T: Analyse de I'expression du genome des levures 
oenologiques en fermentation alcoolique par des approches 
post-genomiques. PhD thesis 2004, Montpellier II University, Sciences des 
Procedes, Sciences des Aliments. 

9. Mendes-Ferreira A, del Olmo M, Garcia-Matinez J, Jimenez-Marti E, 
Mendes-Faia A, Perez-Ortin JE, Leao C: Transcriptional response of 
Saccharomyces cerevisiae to different nitrogen concentrations during 
alcoholic fermentation. Appl Environ Microbiol 2007, 73:3049-3060. 

10. Mendes-Ferreira A, del Olmo M, Garcia-Matinez J, Jimenez-Marti E, Leao C, 
Mendes-Faia A, Perez-Ortin JE: Saccharomyces cerevisiae signature genes 
for predicting nitrogen deficiency during alcoholic fermentation. 
Appl Environ Microbiol 2007, 73:5363-5369. 

11. Contreras A, Garcia V, Salinas F, Urzua U, Ganga MA, Martinez C: 
Identification of genes related to nitrogen uptake in wine strains of 
Saccharomyces cerevisiae. World J Microbiol Biotechnol 201 2, 28(3):1 1 07-1 1 1 3. 

12. Gutierrez A, Beltran G, Warringer J, Guillamon JM: Genetic basis of 
variations in nitrogen source utilization in four wine commercial yeast 
strains. PLoS ONE 201 3, 8(6):e671 66. doi:l 0.1 371/journal.pone.00671 66. 

1 3. Watanabe D, Araki Y, Zhou Y, Maeya N, Akao T, Shimoi H: A loss-of-function 
mutation in the PAS kinase Rim15p is related to defective quiescence entry 
and high fermentation rates of Saccharomyces cerevisiae sake yeast strains. 
Appl Environ Microbiol 201 2, 78(1 1 ):4008-401 6. 

14. Pedruzzi I, Dubouloz F, Cameroni E, Wanke V, Roosen J, Winderickx J, De 
Virgilio C: TOR and PKA signaling pathways converge on the protein 
kinase Riml5 to control entry into GO. Mol Cell 2003, 1 2(6): 1 607-1 61 3. 

15. Swinnen E, Wanke V, Roosen J, Smets B, Dubouloz F, Pedruzzi I, Cameroni E, 
De Virgilio C, Winderickx J: Riml5 and the crossroads of nutrient 
signalling pathways in Saccharomyces cerevisiae. Cell Division 2006, 1(1)3. 

16. Geldermann H: Investigations on inheritance of quantitative characters in 
animals by gene markers II expected effects. TheorAppI Genet 1976, 
47(l):l-4. 

1 7. Sinha H, Nicholson BP, Steinmetz LM, McCusker JH: Complex genetic 
interactions in a quantitative trait locus. PLoS Genet 2006, 2(2):el3. 

18. Sinha H, David L, Pascon RC, Clauder-Munster S, Krishnakumar S, Nguyen M, 
Shi G, Dean J, Davis RW, Oefner TJ, McCusker JH, Steinmetz LM: Sequential 
elimination of major-effect contributors identifies additional quantitative 
trait loci conditioning high-temperature growth in yeast. Genetics 2008, 
180(3;:1661-1670. 

19. Ben-Ari G, Zenvirth D, Sherman A, David L, Klutstein M, Lavi U, Hillel J, 
Simchen G: Four linked genes participate in controlling sporulation 
efficiency in budding yeast. PLoS Genet 2006, 2(1 1 ):el 95. 

20. Deutschbauer AM, Davis RW: Quantitative trait loci mapped to 
single-nucleotide resolution in yeast. Nat Genet 2005, 37(1 2):1 333-1 340. 

21. Gerke JP, Chen CTL, Cohen BA: Natural isolates of Saccharomyces 
cerevisiae display complex genetic variation in sporulation efficiency. 
Genetics 2006, 1 74(2):985-997. 

22. Katou T, Namise M, Kitagaki H, Akao T, Shimoi H: QTL mapping of sake 
brewing characteristics of yeast. J Biosci Bioeng 2009, 107(4)383-393. 

23. Nogami S, Ohya Y, Yvert G: Genetic complexity and quantitative trait loci 
mapping of yeast morphological traits. PLoS Genet 2007, 3(2):e31 . 

24. Kim HS, Fay JO Genetic variation in the cysteine biosynthesis pathway 
causes sensitivity to pharmacological compounds. PNAS 2007, 
104(49):! 9387-1 9391. 

25. Hu XH, Wang MH, Tan T, Li JR, Yang H, Leach L, Zhang RM, Luo ZW: 
Genetic dissection of ethanol tolerance in the budding yeast 
Saccharomyces cerevisiae. Genetics 2007, 1 75(3):1479-1487. 

26. Marullo P, Bely M, Masneuf-Pomarede I, Aigle M, Dubourdieu D: Inheritable 
nature of enological quantitative traits is demonstrated by meiotic 
segregation of industrial wine yeast strains. FEMS Yeast Res 2004, 
4(7):71 1-719. 

27. Smith EN, Kruglyak L: Gene environment interaction in yeast gene 
expression. PLoS Biol 2008, 6(4):e83. 

28. Brauer MJ, Christianson CM, Pai DA, Dunham MJ: Mapping novel traits by 
array-assisted bulk segregant analysis in Saccharomyces cerevisiae. 
Genetics 2006, 1 73(35:1813-1816. 

29. Steyer D, Ambroset C, Brion C, Claudel P, Delobel P, Sanchez I, Erny C, 
Blondin B, Karst F, Legras JL: QTL mapping of the production of wine 
aroma compounds by yeast. BMC Genomics 2012, 13(1):573. 

30. Jara M, Cubillos FA, Garcia V, Salinas F, Aguilera O, Liti G, Martinez C: 
Mapping genetic variants underlying differences in the central nitrogen 
metabolism in fermenter yeasts. PLoS ONE 2014, 9:e86533. 



Page 13 of 14 



31 . Ambroset C, Petit M, Brion C, Sanchez I, Delobel P, Guerin C, Chiapello H, 
Nicolas P, Bigey F, Dequin S, Blondin B: Deciphering the molecular basis of 
wine yeast fermentation traits using a combined genetic and genomic 
approach. G3 2011, 1(4)263-281. 

32. Brion C, Ambroset C, Sanchez I, Legras JL, Blondin B: Variations in 
regulatory networks in yeast revealed in multi-stressed conditions of 
wine fermentation. BMC Genomics 2013, 14:681. 

33. Brem RB, Clinton GYR, Kruglyak L: Genetic dissection of transcriptional 
regulation in budding yeast. Science 2002, 296(5568):752-755. 

34. Flint J, Valdar W, Shifman S, Mott R: Strategies for mapping and cloning 
quantitative trait genes in rodents. Nat Rev Genet 2005, 6(4)271-286. 

35. Keurentjes JJB, Bentsink L, Alonso-Blanco C, Hanhart CJ, Blankestijn-De Vries 
H, Effgen S, Vreugdenhil D, Koornneef M: Development of a near-isogenic 
line population of Arabidopsis thaliana and comparison of mapping 
power with a recombinant inbred line population. Genetics 2007, 
175(2):891-905. 

36. Wang X, Le Roy I, Nicodeme E, Li R, Wagner R, Petros C, Churchill GA, Harris 
S, Darvasi A, Kirilovsky J, Roubertoux PL, Paigen B: Using advanced 
intercross lines for high-resolution mapping of HDL cholesterol 
quantitative trait loci. Genome Res 2003, 13(7):1654-1664. 

37. Lorenz K, Cohen BA: Small- and large-effect quantitative trait locus 
interactions underlie variation in yeast sporulation efficiency. 
Generics 2012, 192(3):1 1 23-1 132. 

38. Satagopan JM, Sen S, Churchill GA: Sequential quantitative trait locus 
mapping in experimental crosses. Star Appl Genet Mol Biol 2007, 
6:Article12. 

39. Fisher RA: The genetical theory of natural selection. In Oxford University: 
The Clarendon Press: 1930. 

40. Lango AH, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, Wilier 
CJ, Jackson AU, Vedantam S, Raychaudhuri S, Ferreira T, Wood AR, Weyant 
RJ, Segre AV, Speliotes EK, Wheeler E, Soranzo N, Park JH, Yang J, 
Gudbjartsson D, Heard-Costa NL, Randall JC, Qi L, Vernon Smith A, Magi R, 
Pastinen T, Liang L, Heid IM, Luan J, Thorleifsson G, ef ah Hundreds of 
variants clustered in genomic loci and biological pathways affect human 
height. Nature 2010, 467(731 7):832-838. 

41 . Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, Cunningham 
MJ, de Andrade M, Feenstra B, Feingold E, Hayes MG, Hill WG, Landi MT, 
Alonso A, Lettre G, Lin P, Ling H, Lowe W, Mathias RA, Melbye M, Pugh E, 
Cornelis MC, Weir BS, Goddard ME, Visscher PM: Genome partitioning of 
genetic variation for complex traits using common SNPs. Nat Genet 201 1, 
43(6)519-525. 

42. Lynch M, Walsh B: Genetics and analysis of quantitative traits. Volume 5. 

In Volume 72. 1st edition. Sunderland, MA01375 USA: Sinauer Associates Inc: 
1998:124. 

43. Manginot C, Roustan JL, Sablayrolles JM: Nitrogen demand of different 
yeast strains during alcoholic fermentation. Importance of the stationary 
phase. Enzyme Microb Technol 1998, 23(7-8)511-517. 

44. Crepin L, Nidelet T, Shanchez I, Dequin S, Camarasa C: Sequential use of 
nitrogen compounds by Saccharomyces cerevisiae during wine 
fermentation: a model based on kinetic and regulation characteristics of 
nitrogen permeases. Appl Environ Microbiol 201 2, 78(22):81 02-81 1 1 . 

45. Hinnebusch AG: Translational regulation of GCN4 and the general amino 
acid control of yeast. Anna Rev Microbiol 2005, 59:407-450. 

46. Ng PC, Henikoff S: SIFT Predicting amino acid changes that affect protein 
function. Nucleic Acids Res 2003, 31:3812-3814. 

47. Daqui T, Li W, Ye Y, Brunger AT: Structure and Function of the 
Yeast U-Box-Containing Ubiquitin Ligase Ufd2p. PNAS 2007, 
104(40):15599-15606. 

48. Benni ML, Neigeborn L: Identification of a new class of negative 
regulators affecting sporulation-specific gene expression in yeast. 
Genetics 1997, 1 47(3):1 35 1 -1 366. 

49. Fendt SM, Oliveira AP, Christen S, Picotti P, Dechant RC, Sauer U: Unraveling 
condition-dependent networks of transcription factors that control 
metabolic pathway activity in yeast. Mol Syst Biol 201 0, 6:432. 

50. Thevelein JM: Signal transduction in yeast. Yeast 1994, 10:1753-1790. 

51 . Donaton MC, Holsbeeks I, Lagatie O, Van Zeebroeck G, Crauwels M, 
Winderickx J, Thevelein JM: The Gapl general amino acid permease acts 
as an amino acid sensor for activation of protein kinase A targets in the 
yeast Saccharomyces cerevisiae. Mol Microbiol 2003, 50(3):91 1-929. 



Brice et al. BMC Genomics 2014, 15:495 
http://www.biomedcentral.com/1471-2164/15/495 



Page 14 of 14 



52. 



53. 



54. 



55. 



56. 



57. 



58. 



59. 



60. 



62. 



63. 



04. 



65. 



66. 



67. 



68. 



69. 



70. 



Grundmann 0, Mosch HU, Braus GH: Repression of GCN4 mRNA 
translation by nitrogen starvation in Saccharomyces cerevisiae. 

J Biol Chem 2001, 276(28):25661 -25671. 

Qui HF, Dubois E, Messenguy F: Dissection of the Afunctional ARGRII 
protein involved in the regulation of arginine anabolic and catabolic 
pathways. Mol Cell Biol 1 991 , 1 1 (4):21 69-21 79. 
Messenguy F, Dubois E: Regulation of arginine metabolism in 
Saccharomyces cerevisiae: a network of specific and pleiotropic proteins 
in response to multiple environmental signals. Food Technol Biotechnol 
2000, 38(4):277-285. 

Phalip V, Kuhn I, Lemoine Y, Jeltsch JM: Characterization of the biotin 
biosynthesis pathway in Saccharomyces cerevisiae and evidence for a 
cluster containing BI05, a novel gene involved in vitamer uptake. 

Gene 1999, 232(1 ):43-5 1 . 

McDonald CM, Wagner M, Dunham MJ, Shin ME, Ahmed NT, Winter E: The 
Ras/cAMP pathway and the CDK-like kinase Ime2 regulate the MAPK 
Smkl and spore morphogenesis in Saccharomyces cerevisiae. Genetics 
2009, 181(2)51 1-523. 

Guldener U, Heck S, Fielder T, Beinhauer J, Hegemann JH: A new efficient 
gene disruption cassette for repeated use in budding yeast. Nucleic Acids 
Res 1996, 24(1 3):251 9-2524. 

Bely M, Sablayrolles JM, Barre P: Automatic detection of assimilable 
nitrogen deficiencies during alcoholic fermentation in oenological 
conditions. J Ferment Bioeng 1990, 70:246-252. 

Sablayrolles JM, Barre P, Grenier P: Design of laboratory automatic system 
for studying alcoholic fermentations in an isothermal oenological 
conditions. Biotechnol Tech 1987, 1:181-184. 
Bezenger MC, Navarro JM: Influence de I'azote sur la fermentation 
alcoolique en milieu modele simulant les conditions de I'oenologie. 
SciAlim 1897, 7:41-60. 

Steinmetz LM, Sinha H, Richards DR, Spiegelman Jl, Oefner PJ, McCusker JH, 
Davis RW: Dissecting the architecture of a quantitative trait locus in 
yeast. Nature 2002, 416:326-330. 

Gresham D, Curry B, Ward A, Gordon DB, Brizuela L, Kruglyak L, Botstein D: 
Optimized detection of sequence variation in heterozygous genomes 
using DNA microarrays with isothermal-melting probes. PNAS 2010, 
107(4):1482-1487. 

Rozen S, Skaletsky HE: Primer3 on the WWW for general users and for 

biologist programmers. In Bioinformatics Methods and Protocols; Methods in 

Molecular Biology. Totowa, NJ: Humana Press; 2000:365-386. 

R Development Core Team: R: A language and environment for statistical 

computing. In R Foundation for Statistical Computing. ; 2012. 

Brem RB, Kruglyak L: The landscape of genetic complexity across 5,700 

gene expression traits in yeast. PNAS 2005, 1 02(5):1 572-1 577. 

Marullo P, Bely (VI, Masneuf-Pomarede I, Pons M, Aigle M, Dubourdieu D: 

Breeding strategies for combining fermentative qualities and reducing 

off-flavor production in a wine yeast model. FEMS Yeast Res 2006, 

6(2):268-79. 

Smyth GK, Speed T: Normalization of cDNA microarray data. Method 2003, 
31:265-273. 

Benjamini Y, Hochberg Y: Controlling the false discovery rate in behavior 
genetics research. Behav Brain Res 2001, 125:289-300. 

Bergstrdm A, Simpson TJ, Salinas F, Parts L, Barre P, Zia A, Nguyen Ba AN, 
Moses AM, Louis EJ, Mustonen V, Warringer J, Durbin R, Liti G: A high- 
definition view functional genetic variation from natural yeast genomes. 

Mol Biol Fvol doi:10.1093/molbev/msu037. 

Liti G, Carter DM, Moses AM, Warringer J, Parts L, James SA, Davey RP, 
Roberts IN, Burt A, Koufopanou V, Tsai LI, Bergman CM, Bensasson D, O'Kelly 
MJ, van Oudenaarden A, Barton DB, Bailes E, Nguyen AN, Jones M, Quail 



72. 



MA, Goodhead I, Sims S, Smith F, Blomberg A, Durbin R, Louis EJ: 
Population genomics of domestic and wild yeasts. Nature 2009, 
458(7236)337-341. 

Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: 
Molecular evolutionary genetics analysis using maximum likelihood, 
evolutionary distance, and maximum parsimony methods. Mol Biol Fvol 
2011, 28:2731-2739. 

Kimura M: A simple method for estimating evolutionary rate of base 
substitutions through comparative studies of nucleotide sequences. 

JMolEvol 1980, 16:111-120. 



doi:1 0.1 1 86/1 471 -21 64-1 5-495 

Cite this article as: Brice ef al:. A genetic approach of wine yeast 
fermentation capacity in nitrogen-starvation reveals the key role of 
nitrogen signaling. BMC Genomics 2014 15:495. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



